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The dynamics of tensile crack fronts restricted to advance in a plane are studied. In an ideal linear 
elastic medium, a propagating mode along the crack front with a velocity slightly less than the 
Rayleigh wave velocity, is found to exist. But the dependence of the effective fracture toughness T(v) 
on the crack velocity is shown to destabilize the crack front if ^ < 0. Short wavelength radiation 
due to weak random heterogeneities leads to this instability at low velocities. The implications of 
these results for the crack dynamics are discussed. 
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The conditions under which solids fail by extension 
of pre-existing cracks are quite well established [Q, but 
the dynamics of crack-front propagation is poorly un- 
derstood. Experiments show a variety of types of crack 
surfaces — including fractal roughness over a wide range 
of length scales [0 — and erratic dynamics of the local 
crack- front velocity ||. With the exception of the uni- 
form motion of a straight crack-front of a planar crack in 
an ideal elastic solid, there has not been, until recently, 
much theoretical analysis of these problems. Indeed, even 
what determines in which direction a section of crack- 
front will advance is not understood |4|-|(|, although it is 
clear that this involves physics on smaller scales than the 
regime of continuum elasticity theory. This has been the 
subject of substantial recent attention |5|||, but the ef- 
fects of heterogeneities, side branching and other effects 
are only just starting to be addressed. 

In this paper, we consider a seemingly much simpler 
problem: the dynamics of a tensile (mode I) crack that 
is confined to a plane in a weakly heterogeneous mate- 
rial. In this case, the local dynamics of the crack-front 
is determined by energy conservation and we can there- 
fore use continuum mechanics with the non-linear small 
scale physics near the crack-front just determining the 
local fracture toughness. But we shall see that, never- 
theless, the combined effect of even weak heterogeneities 
and elastodynamics can be dramatic. 

In a scalar approximation to elasticity, Perrin and Rice 

found that to lowest order in a weak randomness per- 
turbation theory, heterogeneities in the fracture tough- 
ness led to slow — logarithmic in time — roughening of a 
planar crack-front suggesting that this system does not 
reach a statistical steady state. In this paper, we first 
show that the dynamics of a tensile crack are very differ- 
ent from the scalar approximation. In particular, in an 
ideal system with constant fracture toughness, T, a mode 

1 crack that is advancing with velocity v, supports a prop- 
agating mode along the crack-front with a velocity that 
is slightly slower than the Rayleigh wave velocity, c, and 
depends weakly on v. This mode does not exist in the 
scalar approximation and its presence dramatically en- 
hances the effects of weak heterogeneities. Furthermore 
if the toughness T(v) is not constant but decreases with 



velocity, the propagating mode — and hence the crack- 
front — becomes unstable. A perturbative analysis of the 
effects of random local heterogeneities is found to yield 
an effective T(v) whose contributions from short wave- 
length acoustic emission indeed make ^ < at least 
for v/c small. The implications of these results are dis- 
cussed and we speculate on possible scenarios suggested 
by them. 

We study a tensile crack confined to the plane y = 
with the crack open in the region x < F(z,t). The nor- 
mal displacement has a discontinuity across the crack, 
i.e., u y (y = + ) ^ u y (y = 0~), while the normal stresses, 
a iy(y = C^) vanish on the crack surfaces. Under an ex- 



ternal load, the displacement field u y has a y/F{z,t) — x 
singularity at the crack-front with an amplitude propor- 
tional to the local mode I stress intensity factor, Ki(z, t). 
As long as the crack remains planar, the loading is purely 
mode I so we write K = Kj. The system is under a static 
load applied very far away so that for a straight crack 
at rest, (i.e., F=constant), K — X°°=constant. The 
crack is considered to be initially at rest (at large nega- 
tive times) and move a total distance small compared to 
its length and the sample dimensions and for a time less 
than the sound travel times across the sample. 

If the crack-front advances, F — > F+SF, an energy per 
unit area of new crack T[x = F{z, t), z] must be provided 
to the crack front in order to fracture the solid; in an ideal 
quasi-equilibrium situation this is just twice the solid- 
vacuum inter-facial energy density, more generally it is 
the local fracture toughness that includes the effects of 
small scale physics for which linear continuum elasticity 
is not valid. The fracture energy will be provided to the 
crack- front by a flux of stored elastic energy per unit area 
of new crack G{z, t, {F}), which in general depends on 
the past history of the crack and its instantaneous local 



velocity ^ . The dynamics of the crack is given by 

G[z, t, {F{t' < t)}} = T[x = F(z, t), z] 

for all z and t. 

The available energy, G, has the general form 

G = $[v ± (z,t)] G R [z,t,{F(t' <*)}] 



(1) 



(2) 
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where Gr, which is independent of dF(z, t)/dt), is the 
available elastic energy at (z, t) for a crack at rest locally 



and ^-L = ^|f"/yl+(fjf") is the local velocity normal 
to the crack front. For a straight crack F =constant, 
G R = G°° = (K°° f (1 - u 2 ) JE, with E Young's mod- 
ulus and v Poisson's ratio || . The fraction of Gr avail- 
able for fracture decreases linearly from unity 
for small v± and goes to zero for vj_ = c, yielding, for 
a straight crack with uniform toughness, a monotonic 
v{G°°) = tp-^G^/T) for G°° greater than the Griffith 
threshold, i.e., G°° > T. 

In a weakly inhomogeneous material, one can hope to 
perturbatively expand about a uniformly opening crack, 
writing 



F(z,t)=vt + f{z,t), 
T(x,z) = T [l + j(x,z)] 
and G = G [l + g(z, t)]. 



(3) 



The variation in the available energy for a small crack- 
front distortion has the form: 

g = -P®f + S(f, /) + T(f, /, /) + 0(/ 4 ) (4) 

where in Fourier space z — > k,t — » lu 

2c 



P=\k\{ -T— J^-T 2 -- 
&■ — v z a 

+ — I dWI(W,a 2 ) } 
2m J 



(5) 



with 



2v 2 W -a 2 (W + v 2 ) 
^/W{W - a 2 )(W - v 2 ) 2 

2 



xln 



b 2 



b 2 



(6) 



the contour integral circling the cut in the complex W 
plane from W = b 2 to W — a 2 with a and b the longi- 
tudinal and transverse sound velocities respectively; the 
argument of the logarithm in Eq. (^) the function whose 

zero is W = c 2 ; a 2 = w + v 2 the square of a phase ve- 
locity in the rest frame; and lu — > lu + iO needed to define 



all cuts e.g., Vc 2 — u 2 



-isign(u})^/ a 2 



c for large lu. 



The functionals S and T are bi- and tri- linear, respec- 
tively, with each consisting of a sum of many terms that 
are sequences of multiplications and convolutions [®) in- 
volving P(k,uj;v) and functions related to it The linear 
kernel, P, can be derived from the results of Movchan and 
Willis H but a simpler method was developed to obtain 
the non- linear terms; it will not be reproduced here. For 

large uj, P w — iuj (— ^ representing the dependence 

of G on local velocity, while for w — 0, P oc |A;| repre- 
senting the long range stiffness of the crack (as found in 
the scalar elasticity approximation) [Q. Generally, for 



\ui\ < v 'c 2 — v 2 \k\, the imaginary part P"(k,ui) — since 
no bulk or crack-surface modes can carry away energy in 
this regime. But, perhaps surprisingly, for a particular 
lu < c\k\, P[uj — ±s(v)\k\] = implying a propagating 
mode along the crack-front with velocity s(v) in the mov- 
ing frame and , s i v ^ * increasing from 0.94 for v — to 1 
for v — > c. Evidence for this propagating mode is seen in 
recent numerical simulations of Morrissey and Rice jio) . 
This is in striking contrast to the scalar approximation 
for which P = c i°i v '± \Jk 2 {c 2 — v 2 ) — uj 2 yielding a 1/Vt 
decay of the maximum distortion at long times after a 
localized disturbance caused by a tough region. 

Damping of bulk sound modes can be taken into ac- 
count by replacing uj 2 by lu 2 / (1 — iLUTd) in P (and S and 
T) . This yields a damping rate of the propagating crack- 
front mode proportional to tu 2 Td- The leading effect of 
weak heterogeneities can readily be seen to yield a mean 
square roughness in steady state: 

C(z)^(\f(z + z',t)-f(z',t)f) 

f (q= -^,fc)(2-2cosfcz) 
v\P{k,w)\ 2 



(7) 



with T the Fourier transform of the correlations 
(7(2;, z) r y(x', z')) of the fractional toughness variations. 
The propagating mode yields at large separations \z\^> 
cl/v 



C(z) 



M 
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(8) 



where I is a characteristic correlation length of the tough- 
ness variations. This is much rougher than the C(z) ~ 

^ In (\z\/cTd) or C{z) ~ ffi Qz\/cr,i) obtained, re- 
spectively, in the scalar or quasi-static (P = —iw/c+ \k\) 
approximations to the stress transfer. 

As we shall see, the roughness of Eq. (B) is caused 
by a straight moving crack being at the boundary of a 
regime of stability so that small perturbations — which 
we will find can be driven by the randomness induced 
fluctuations — can cause it to become much more unsta- 
ble. 

We are particularly interested in the behavior on 
length and time scales much longer than the scales of 
variations of the toughness, i.e., I and l/v respectively. 
On long scales, one might expect that the crack-front 
would have a well defined coarse-grained position / that 
does not vary much, even if on small scales the local ve- 
locities and local angles, djjdz, are strongly fluctuating. 
The long wavelength and low frequency components of 
the elastic energy flux into the crack, G (or g) will still 
be given by g = —P <S> f ■ But the non-linear effects of 
the random toughness and the non-linearities in G[{F}] 
will make the effective toughness T include contributions 
from radiation of short wavelength sound waves as well 
as other small scale physics such as plasticity, etc. These 
cause the effective toughness to increase and become de- 
pendent on the coarse-grained velocity and history. If the 
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randomness is effectively averaged out, on long scales the 
crack will be moving almost uniformly with a mean ve- 
locity v determined by 



f[v, {r}] = g°°^(v) 



(9) 



with r depending on the distribution and spatial corre- 
lations of the microscopic T(x, z). 

Naively, it would appear that this system is stable as 
long as v{G°°) is single valued, i.e., if A 



> 



dv 



But, in fact, the elastodynamics implies that it is unsta- 
ble at finite wavelengths unless ^= > 0. Linearizing in 

perturbations / about a putative fixed velocity steady 
state, given by Eq (||), we have for long wavelengths and 
low frequencies, —iujAf = —Pf. For A > 0, the response 



function \ 



has no singularities in the upper half 



plane and thus the system is stable. The perturbative ef- 
fects of the residual small randomness will then result in 
fluctuations in / that are strongly suppressed and more 
like the quasi-static result, with C{z) ~ -~ lnz for small 
A. 

With the dynamics appropriate with A > 0, it can be 
shown that non-linear terms consistent with the symme- 
try are irrelevant at long scales (in contrast to the case 
of a simple elastic line moving through a random poten- 
tial). Thus the long wavelength low frequency behavior 
for A > is quite simple, although the dynamic correla- 
tions will still reflect the elastodynamics unless A is very 
large. But for A < 0, \ will have poles in the upper half 
plane corresponding to unstable propagating modes with 
u> S3 ±s(u)|fc| — irj{v)\k\ with rj oc Ac 2 < 0. 

Since we have begun with a more microscopic descrip- 
tion of the system, it should be possible to compute the 
large scale response function perturbatively for weak het- 
erogeneities and thereby find the sign of A for the case 
in which the microscopic T(x, z) is not velocity depen- 
dent. It is somewhat simpler to compute a slightly differ- 
ent quantity: x{k,u>) = S /e(k'u) wn ich is the response to 
extra energy per unit area added at the crack-front (or 
roughly equivalently a decrease in the toughness); this 
will yield the behavior of perturbations and hence also 
the randomness induced fluctuations. 

The unperturbed response function is x — Xo = -p- 
Since the interesting regime is near where P vanishes, a 
direct expansion of x m powers of T is not useful. Thus 
instead, we expand 



X 



(10) 



with E (the analog of a "self-energy" ) expected to be well 
behaved at small k and u>; we will thus focus on this. 
Conventional Feynman diagrammatic techniques can be 
used via a Martin-Siggia-Rose formulation of the dynam- 
ics, analogous to treatments of, e.g., charge density wave 
dynamics jllj]. The main complication here is the form 
of the non-linearities. In addition to the non-linearities 
from expanding T[x = vt + f(z, t)] in powers of /, there 



are also the non- linear terms in G[{F}], Eq. (|J), aris- 
ing from the moving boundary problem of an irregular 
crack- front. To compute even the leading correction to 
X -1 , both the second and third order non-linearities in 
G are needed. But for velocities v <C c, only the second 
order parts, S, are needed and there are some simplifica- 
tions; we will focus on this regime here. 

From the structure of the perturbative expansion, it 
can be seen that E(fc, lu) vanishes for k — u> — 0, and for 
small fc and lj has the form 



E(fc,w) S3 iuA + BP(k,w) 



0(Lu 2 ,k 2 ,P 2 ,P 2 ,tLU 3 ...) 

(11) 



with P2 = J" <ko' 9 ' . From the above discussion, it 
should be clear that the coefficient A will be the most 
important, with A S3 A[1+0(B)]. We are thus interested 
primarily in the weak randomness expansion of A. 

For low velocities, inverse powers of v are associated 
with non-linearities from the random toughness and thus 
would be expected to dominate, with the natural expan- 
sion parameter being {^ 2 )c/v But the leading low 
velocity contribution to the expansion of A is found to 
vanish exactly. The dominant term thus involves both 
the non-linearities from the random toughness and those 
from the changes in crack shape. Evaluating these terms 
yields 



A: 



dq\q\ 
2tt 



T(g,fc = 0) 



duj\k\ [ImP(k,u)Y 



2ttw 2 



\P(k,u) 



(12) 



which is independent of v. Corrections to this arise pri- 
marily from the randomness induced non-linearities with 
(j 2 ) 2 /v 2 terms expected at next order. Thus the expan- 
sion should be good if (j 2 )<^iv 2 /c 2 (rather than the naive 
( 7 2 )<i;/c). 

At low velocities, we thus conclude that the crack- 
front is unstable at long wavelengths, A > A c . The 
iu 3 term in and the sound wave damping stabi- 
lize the system for shorter wavelengths, A < A c , with 



A, 



CT d 

W) 



'(f) 



3/2 



note that the second cutoff 



corresponds to modes with period much longer than the 
time, -, for the crack to advance a distance I, The most 
rapidly growing modes will have A ~ 2A C . 

At higher velocities, there are contributions of both 
signs to the order 7 2 terms in E and the sign of A de- 
pends on more details of the correlations, Y(p, k), but 
at least for simple forms of the correlations, A can be- 
come positive for high velocities and the crack-front will 
be stable. @ 

We now consider various possible scenarios for mode 
I planar cracks. The simplest scenario, suggested by 
the perturbative analysis, is that there will be a regime 
of unstable velocities probably extending down to zero 
velocity — although the perturbative analysis certainly 
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fails there. Associated with this, one might expect hys- 
teresis in the macroscopic velocity v versus the load K°° 
(or G°°) with a stationary crack jumping to a finite ve- 
locity as the load is increased, with a statistically steady 
state moving phase that is only logarithmically rough — 
corresponding to A(v) > 0. On decreasing the load, the 
crack would not stop until it reaches some lower criti- 
cal load at which A — > whereupon it would stop sud- 
denly. Alternatively, for macroscopic systems, the hys- 
teresis might disappear leaving a single-valued but dis- 
continuous v(G oa ) and a strongly fluctuating crack for 
G°° just above the critical load Gf « T[l + 0((j 2 })}. 
The scenario with the stationary crack jumping to a fi- 
nite velocity when the loading is increased, is consistent 
with our numerical results on a simplified model, fill 

Another possibility is that strongly non-linear fluctu- 
ations in the local velocity driven by the variable tough- 
ness might drive A positive and stabilize the crack dy- 
namics even for low velocities. But our understanding 
of the behavior near G£° in the quasi-static approxi- 
mation suggests that this is unlikely |14|. Perhaps the 
most interesting possibility is that for some range of ve- 
locities, a non-trivial, non-linear statistical steady state 
might occur that is much rougher than logarithmic and 
perhaps characterized by propagating modes along the 
crack. This would be roughly analogous to the behavior 
of the KPZ fl5| equation in dimension d > 2 for which 
both a non-trivial strongly fluctuating phase and a simple 
weakly fluctuating phase exist. Combinations of these 
possibilities could also occur. The behavior should de- 
pend on aspects of the correlations of the random tough- 
ness as well as its magnitude and on velocity dependent 
toughness that arises from other physics such as plastic- 
ity, dislocation emission, thermal creep, and small scale 
side branching. 

Note that even in an ideal system with no tough- 
ness variations, thermal fluctuations can give rise to ef- 
fects similar to those discussed here as the elastic energy 
needed to break a bond is less if the bond is already 
stretched due to thermal fluctuations. For v a substan- 
tial fraction of c, the growth rate of the fluctuation driven 
instability will be of order i<5 4 with S of order the ratio 
of the temperature to the bond breaking energies, and r 
a microscopic time. In practice, the effects of the atomic 
scale physics on the effective T(v) will probably dominate 
over thermal fluctuation effects. 

One type of correlation in the toughness that is impor- 
tant to consider is periodic inhomogeneities — a crude ap- 
proximation to some of the effects of a crystal lattice. In 
this case, weakly variable toughness can again be shown 
to lead to a perturbatively unstable crack-front at low 
velocities. But the short wavelength acoustic radiation 
that leads to the velocity dependent effective toughness 
arises from high order harmonics — of order = — of the 
periodicity in the x-direction; thus at low velocities the 
leading perturbative effect can be strongly suppressed. 
This instability caused by non-steady motion in a peri- 
odic system would appear to be related to that found by 



Marder and Gross in a one-dimensional "ball and spring" 
model of a crack Jll| . 

We have shown here that both random and periodic 
inhomogeneities in the toughness can dynamically desta- 
bilize the front of a tensile crack moving in a plane and in 
general such cracks will be unstable to finite wavelength 
disturbances if their effective toughness decreases with 
velocity. Understanding the resulting dynamic states of 
the system as well as the more challenging problem of 
out-of-plane crack motion, we must leave for future work. 

We would like to thank J. R. Rice and J. Morrissey 
for useful discussions and communicating their numer- 
ical results. This work has been supported in part by 
the NSF via DMR-9106237, 9630064 and Harvard Uni- 
versity's MRSEC. 



[1] B. Lawn, Fracture of Brittle Solids, Second Edition, 
Cambridge University Press, Cambridge 1993. 

[2] P. Daguier, S. Henaux, E. Bouchaud and F. Creuzet, 
Phys. Rev. E 53, 5637 (1996) and references therein. 

[3] J. Fineberg et al, Phys.Rev.Lett. 67, 457 (1991). 

[4] An account of various proposed criteria is given in 
H. Bergkvist and L. Guex, Int. J. Fracture 15, 429 (1979). 

[5] J. Hodgdon and J. P. Sethna, Phys. Rev. B 47, 4831 

(1993) ; B. Cotterell and J. R. Rice, Int. J. Fracture 16, 
155 (1980). 

[6] E. S. C.Ching, J. S. Langer and H. Nakanishi, Phys. Rev. 

Lett. 76, 1087 (1996). 
[7] G. Perrin and J. R. Rice, J. Mech. Phys. Solids, 42, 1047, 

(1994) . 

[8] L. B. Freund, Dynamic Fracture Mechanics, Cambridge 

Uni versity Press, Cambridge, 1990. 
[9] J. R. Willis and A. B. Movchan, J. Mech. Phys. Solids 

43, 319 (1995). 

[10] J. Morrissey and J. R. Rice, Trans. Amer. Geophys. 
Union 77, F485, (1996). 

[11] By analogy with charge density waves; see Onuttom 
Narayan and Daniel S. Fisher, Phys.Rev.B. 46, 11520 
(1992), and references therein. 

[12] We note that in the scalar model rather different physics 
occurs: the short wavelength radiations lead to < 
resulting in a supersonic mode at low velocities . This 
mode is generically weakly unstable, resulting again in 
the instability of the moving crack-front. At higher veloc- 
ities, depending on the microscopic details of the system, 
can change sign resulting in stable crack-fronts that 
are only logarithmically rough 

[13] S. Ramanathan and D. S. Fisher, Unpublished. 

[14] If A is large and positive , fluctuations will drive it to- 
wards low velocities, as can be shown following Ertas 
and Kardar, Phys. Rev. E 49, R2532 (1994) and Ref. 
jllj|. This is likely to make the crack-front unstable. 

[15] M. Kardar, G. Parisi and Y-C. Zhang, Phys.Rev.Lett, 
56, 889 (1986). 

[16] M. Marder and S. Gross, J. Mech. Phys. Solids, 43, 1, 

(1995) . 



4 



